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I. INTRODUCTION 

Understanding the origin of ferromagnetism was one 
of the success stories of applying quantum mechanics to 
solid-state systems. Classically, magnetic moments inter- 
act via the dipole-dipole interaction, which is much too 
weak to explain magnetic order at as high a temperature 
as is observed in, e.g., iron or nickel.^ The explanation 
of this conundrum was found to be the exchange inter- 
action mechanism, which leads to a spin-spin interaction 
that is governed by the Coulomb interaction via the Pauli 
principle. This was first understood in the context of 
atomic and molecular physics in the 1920s, and applied 
to solid-state physics in the 1950s. ^ Somewhat ironically, 
a straightforward application of the exchange interac- 
tion concept leads to a spin-spin interaction that is too 
strong, as the relevant energy scale is the atomic scale, 
or roughly 100, 000 K. Many-body and band-structure 
effects renormalize this scale and bring it down to the 
observed ferromagnetic scale of rougly 1, 000 K or lower 
This is still much larger than the dipole-dipole scale, and 
the latter is often neglected in the discussion of ferromag- 
nets. When it is considered, e.g., for its influence on the 
critical behavior,^"— it is usually added phenomenologi- 
cally to models that describe the exchange interaction. 
Another spin-spin interaction that has been of interest 
lately is the Dzyaloshinsky-Moriya (DM) interactioui^ 
It results (in systems with suitable lattice structures) 
from the spin-orbit interaction, has been derived from 
microscopic models, and is believed to be responsible 
for the helical magnetic order observed in MnSi and 
FeCe.® Rough estimates show that the DM interaction 
and the dipole-dipole interaction are of about the same 
strength, and should thus be considered togetheri^ Fur- 
thermore, in weak ferromagnets, which order only at low 
temperatures, all three interactions can be comparable 
in strength, which can make the dipole-dipole and DM 
interactions crucial. 



In this paper we provide a comprehensive derivation 
of all of these effects within one unified framework, 
namely, a field-theoretic description of electrons and pho- 
tons. Starting with finite-temperature quantum elec- 
trodynamics (QED) coupled to a field-theoretic descrip- 
tion of finite-density quasi-relativistic electrons we show 
that the exchange, dipole-dipole, and DM interactions 
all appear naturally upon integrating out the photons. 
The exchange and DM interactions arise from integrating 
out the scalar part of the electromagnetic potential; the 
dipole-dipole interaction, from integrating out the vector 
potential. Furthermore, the DM and dipole-dipole inter- 
actions are indeed of the same order in the relativistic 
corrections to the Schrodinger equation (i.e., of second 
order in v-p/c, with v-p the Fermi velocity and c the speed 
of light, or of second order in the fine structure constant 



Integrating out the fermions then leads to an ef- 
fective theory for quantum magnets that general- 
izes and replaces the Hertz-Millis theor y^°i^^ and its 
generalizations.™ More generally, the theory provides a 
derivation of spin-spin interactions in itinerant Fermi sys- 
tems in general, whether or not they are in a parame- 
ter regime where they develop long-range magnetic or- 
der. Our results are therefore relevant, for instance, for 
fermionic atoms in optical traps or on optical latticesJ^ 

This paper is organized as follows. In Sec. we con- 
sider, as a warm-up and to introduce various concepts, 
classical magnets, and show how the vector potential cou- 
pling to the magnetization gives rise to the dipole-dipole 
interaction. In Sec. IIIII we develop the technical machin- 
ery for dealing with quantum magnets and provide the 
derivations mentioned above. In Sec. IIVI we discuss our 
results and provide a summary and conclusion. Some 
technical details are relegated to various appendices. 
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II. EFFECTIVE ACTION FOR CLASSICAL 
FERROMAGNETS AND HELIMAGNETS 

We now proceed to derive an effective action for mag- 
nets that includes the effects of the fluctuating electro- 
magnetic potential. We first consider the classical case 
as a warm-up; we will generalize to the quantum case in 

secnni 



A. Dipole-dipole interaction 

Consider a classical model for a ferromagnet with a 
three-component order parameter M. In addition to the 
field M{x) we need to consider the electromagnetic vec- 
tor potential A{x), and the partition function Z is given 

by 



Z = J D[M,A] e^I^^^l 
'd[M] e--^[^l/^. 



(2.1a) 
(2.1b) 



The model is defined by specifying the action 5, and in 
Eq. (j2.1bl) we have anticipated integrating out the vec- 
tor potential to obtain an effective action T in terms 
of the order parameter only. T denotes the tempera- 
ture, so J- is the free energy in mean-field approxima- 
tion. Throughout this paper we will use units such that 
Boltzmann's constant and Planck's constant are equal to 
unity, kB ~ h = I. 

For the order-parameter part of 5, we consider an 
0(3)-symmetric (/)'*-theory, 



Sm - -jr 



dx 



iM^x) + ^{VM{x)r + ^M\x) 



(2.2a) 

Sm represents a Landau-Ginzburg- Wilson (LGW) theory 
of an isotropic ferromagnet with volume — >■ cxd. The 
parameter t contains the exchange interaction that leads 
to a magnetic ordering transition. In mean-field approxi- 
mation this transition occurs at i = 0, with t > describ- 
ing the paramagnetic phase, and t < the ferromagnet- 
ically ordered one. a > and m > are two additional 
model parameters, and (ViW)^ = di Mj d"^ . Here, 
and throughout the paper, summation over repeated vec- 
tor, tensor, and spinor indices is implied unless otherwise 
noted. Note that Sm is separately invariant under rota- 
tions in M (spin) space and real space, respectively. 

The magnetizationii M couples linearly to the curl of 
the magnetic vector potential A: 



[ dx M{x) ■ (V X A{x)), (2.2b) 
Jv 



with /iB — e/2meC the Bohr magneton in terms of the 
electron charge e, the electron mass rrie, and the speed of 
light c. A and V x ^ transform as vectors in real space, 
and therefore Sc is invariant only under co-rotations of 



spin space and real space. It is this coupling of the mag- 
netization to the fluctuating vector potential that allows 
one to consider the magnetization as having a particular 
direction in real space. The vector potential is governed 
by 



Sa = 



-1 



SttT 



dx 



(V X A{x)f 



A{x)f 



(2.2c) 



with p any real number. The first term in Eq. ()2.2c|) is 
the magnetic energy, and the second term with coupling 
constant 1/p is a gauge fixing term. One popular choice 
is p = 0, which enforces a Coulomb gauge, V • A = 0; 
another one is the Feynman gauge, p = Iji^ii^ Either 
choice ensures a finite A-propagator. In Coulomb gauge, 
it is 



{A,{k)A,{-k)) 



fc2 



(2.3) 



The vector potential can now be integrated out exactly, 
which leads to an effective action in terms of M only. 
Alternatively, we can consider the magnetic induction 
B — 'V xA the fundamental field to be integrated out. In 
that case, the gauge fixing condition needs to be replaced 
by a constraint that enforces the Maxwell equation V ■ 
B = 0. That is, the Eqs. (g^B, (l2:2cl) . and (EH) are 
replaced byii 



Sa 



SttT 



T 



dx 



/ dx M{x) ■ B{x), 
Jv 



B^{x) + -{V-B{x))^ 



{B,{k)B,{-k)) = A^T 



(2.2b') 



p-s-O 

(2.2c') 



(2.3') 



Either way we find 



J" = 



dx 



--M\x) 
4 ^ ' 



\{t-A^pl) M^{x) + ^{VM{x)f 
+ 27r/i| dij{k) Mi{k) Mj{-k), (2.4a) 



with 



dij (fc) — kikj 



(2.4b) 



The terms generated by integrating out the vector poten- 
tial we recognize as the leading contribution to the dipolc- 
dipole interactional^ plus a shift of the Landau parameter 
t by 47r/i|. The scalar potential ip{x), whose gradient is 
the electric field, does not lead to any magnetic interac- 
tions in a classical theory. This changes once the system 
is treated quantum mechanically, see Sec. IIIII below. 
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FIG. 1: A diagram that generates a (V ■ M)'^ term. 



B. Renormalization, and higher order terms 



The dipole-dipole operator dij, Eq. (j2.4bp . transforms 
as a rank-two tensor in momentum (or real) space, 
and M transforms as a vector in spin space. Conse- 
quently, the dipole-dipole interaction is invariant under 
co-rotations in real space and spin space. This raises the 
question of other terms in the action that have the same 
symmetry properties. For instance, (V ■ M)'^ is allowed 
by symmetry. This term, and terms of higher order in 
the gradient, are generated by a renormalization of the 
action T, as we now proceed to show. 

A renormalization of the action J- generates additional 
terms by virtue of the anisotropic M-propagator, which 
now reads 



{A'Uik) M,{-k)) 



(2.5) 



For instance, two-loop diagrams of the structure shown 
in Fig.[T]both renormalize dij (k) and lead to a new vertex 



J2hkjM,{k)Mj{-k), 



(2.6) 



as well as to higher order anisotropic gradient terms. 
Equation (|2.6I) represents the (V • M)'^ term that was 
mentioned above. Conversely, if one starts with a theory 
that contains a (V • iW)^ term, which is allowed by sym- 
metry and hence should be included in any Landau the- 
ory, then the leading dipole-dipole term will be generated 
in perturbation theory even if it was not included in the 
bare action. The complete LGW action for a classical, 
isotropic Heisenberg ferromagnet, up to terms quadratic 
in gradients and quartic in the order parameter, thus 
reads 



T = 



dx 



-M'{x) + -{VMix)r + -M\x) 



d f 

- y J dxdy M, (x) d,j {x ~ y) Af, (y ) 

+ ^ / dx {V-M{x)f . (2.7a) 

Here dij{x — y) is the Fourier transform of dij{k) in Eq. 
(|2.4bp . namely, 

dij (x-y) = 



1 



dxi dxj \x — y\ 



(2.7b) 



r is the bare distance from the critical point in the effec- 
tive LGW theory that takes into account the effect of the 



vector potential, and a, u, do, and d2 are the remaining 
Landau parameters. As is clear from the above discus- 
sion, one expects the bare values of do and c?2 to be small 
of order 1/c^ compared to the other parameters in nat- 
ural units. These two parameters are usually set equal 
to zero in elementary treatments of classical Heisenberg 
ferroniagnets. 



C. Dzyaloshinsky-Moriya interaction, and 
helimagnets 

The terms in the action so far are all even in the gradi- 
ent operator, and hence invariant under spatial inversion. 
The spin-orbit interaction can eliminate this requirement 
by coupling the electron spins to the underlying lattice, 
provided the crystal structure is not inversion invariant. 
Dzyaloshinsky^ and Moriya^ (DM) showed that to linear 
order in the spin-orbit interaction gso the relevant term 
is 



J'DM — 



-cdm 
2T 



dx M{x) ■ (V X M{x)) , 



(2.8) 



with cdm oc gso- At a classical level, this term is purely 
phenomenological. DM showed how to derive it in the 
context of quantum mechanics, and in Sec. IHII we will 
see how it arises in the context of field theory. 



III. QUANTUM SYSTEMS 

We now turn to a quantum mechanical description of 
itinerant fermion systems in general, and certain types of 
magnetism in particular. We will show how the exchange 
interaction, the dipole-dipole interaction, and the DM in- 
teraction naturally arise in the context of a field-theoretic 
description of itinerant electrons. The former two lead to 
ferromagnetism, and the latter, if it is present, to heli- 
magnetism. 



A. Action 

In Appendix |X] we list the complete action for free 
quasi- relativistic electrons, to order coupled to the 
electromagnetic field. Several terms in the complete ac- 
tion are not relevant for our present purposes. The only 
effect of the Darwin term, Eq. ( |A4g[ ), is to modify the 
Coulomb interaction on length scales given by the elec- 
tronic Compton wave length Ao = 1/tocC. The relativis- 
tic mass enhancement, the last term in Eq. (IA4bl) . is 
a higher-order gradient term that is small compared to 
terms of the same form that are generated by renormaliz- 
ing the final effective action. Finally, the Landau diamag- 
netic terms, Eq. (j A4f P . give rise to diamagnetism and, in 
the presence of an external magnetic field. Landau levels. 
These effects are physically very different from ferromag- 
netism or helimagnetism, and we do not consider them 
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here. Finally, in quantum electrodynamics the Fadeev- 
Popov ghost field does not couple to any other fields. Its 
only effect is to subtract the contribution of the unphys- 
ical longitudinal photon polarization to the free energy, 
and one has to keep it only if one is interested in the ab- 
solute value of the latter. Neglecting all of these terms, 
we thus consider the following action: 



S[ip,ip;Af_,] ^ dx tpa{x) 



-dr 



1 



2mo 



dx f(x) n(x) 



^J■B J dx B[x) 



■ ns{x) 



4y„2^2 / dx 1pcri{x) Crcria2 



(Vifix) X V)^P,,ix). 



(3.1) 



Here x = {x,t) comprises real space position x and 

imaginary time r, and J dx = Jy dx J^^^ dr. and 
are Grassmann-valued fields for electrons with spin 
projection cr, and the first term in Eq. ()3.ip describes 
free electrons with chemical potential ji}^ = (0, —Ai) 
(/J, = 0, 1, 2, 3; j = 1, 2, 3) denotes the electromagnetic po- 
tential, with (f) the scalar potential and A = (Ai, A2, A3) 
the vector potential, and the second term in in Eq. (13. ip 
describes the free electromagnetic field in Feynman gauge 
(i.e., p = 1 in Eq. \AAc\i )^ Note that both 4- vector po- 
tential fields in the second term carry covariant indices; 
that is, the A-action is euclidian.— n{x) = ipa{x) ipa{x) 
and rig (a;) — ipa-i {x) Ua^a^ ipa-^ (x) are the electronic num- 
ber and spin density^^ respectively, with a = {a^^a^, a^) 
the Pauli matrices, and the third term in Eq. de- 
scribes the coupling of the electrons to the electromag- 
netic field, with S = V X A the magnetic induction. 
Finally, the last term in Eq. (|3.ip describes the spin-orbit 
interaction. Note that both terms coupling the scalar po- 
tential Lp to the ferniions carry an extra factor of i com- 
pared to what one might expect from the first quantized 
Hamiltonian. This has the same origin as the Euclidian 
metric mentioned above.=22. 

Equation p.ip describes a continuum model. Some 
of the effects we are interested in are present only in the 
presence of certain types of lattices, and we will comment 
later on the modifications that occur if the electrons are 
put on a lattice. 



B. Integrating out the photons 



The action, Eq. ()3.1|) . depends only bilinearly on the 
electromagnetic potential. The latter can therefore be 
integrated out exactly, albeit at the expense of creat- 
ing four-fermion terms. The latter represent electron- 
electron interactions that are mediated by the exchange 
of virtual photons. Technically, we need the photon prop- 
agator, which we can read off the second term in Eq. 



{A^{x)A,{y))=6^,,V{x-y), 



(3.2a) 



with 



p(a:-y) = -47r (a^/c^ + V^) ^6{x-y), (3.2b) 



or, in Fourier space, 

P(fc)=47r/ {nl/c^ + e) 



(3.2c) 



Here k = (ir2„,fc) comprises a bosonic Matsubara fre- 
quency ri„ = iirTn and a wave vector k. 

The result of integrating out the photons exactly 
is very complicated and involves interacting electronic 
modes in both the spin-singlet and spin-triplet channels, 
the particle-particle and particle-hole channels, and all 
angular-momentum channels. We will restrict ourselves 
to those terms that are most relevant for magnetism, i.e. 
interactions between spin-density fluctuations, or modes 
in the s-wave particle-hole spin-triplet channel. 



1. 0{1/(P): Exchange interaction 

We organize the various contributions to the effective 
electron-electron interaction in powers of 1/c. To zeroth 
order only the term coupling ip to the number density 
n in Eq. p.ip contributes. Integrating out Lp leads to a 
Coulomb interaction 



with 



i / dxdy n{x)vc{x-y)6{Tx 



vc{x) = e^/\x\. 



)n{y), (3.3a) 



(3.3b) 



Here we have neglected the dynamical nature of the pho- 
ton propagator V and have replaced it by its value at 
ri„ = 0. The reason for this approximation is that 
Fermi-liquid effects lead to a dynamical screening of the 
Coulomb interaction that is a much larger effect than the 
relativistic dynamics inherent in Eqs. p.2|) . 

Equation (|3.3ap contains number density fluctuations 
at all wavelengths. If one restricts the theory to inter- 
actions between long-wavelength fluctuations, then this 
interaction can be rewritten as a sum of parts that in- 
cludes an interaction between spin-density fluctuations, 
see Ref. [l^ and Appendix [BJ The basic point is that an 
interaction between number density fluctuations at large 
wave numbers can be written as one between spin den- 
sity fluctuations at small wave numbers. In an effective 
low-energy theory that contains only fluctuations at wave 
numbers smaller than some cutoff A, Eq. p.3ap therefore 
contains a contribution 

Sex = Y / dx ns{x) ■ ns{x), (3.4) 
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where the prime on the integral indicates that only the 
small- wave-number contributions (smaller than A) to the 
spin density rig are to be considered in order to avoid 
overcounting. As has been explained in Ref. '23', it is 
convenient to choose the cutoff A as a fixed fraction of 
the Thomas-Fermi screening wave number, and the spin- 
triplet interaction amplitude Ft is a Fermi-surface aver- 
age over wc(fc — p) Q{\k — p\ — A), with k and p pinned 
to the Fermi surface. The restriction to small wave num- 
bers will be understood from now on, and we will drop 
the prime on integrals. 

5'ex is the exchange interaction between electronic spin- 
density fiuctuations that leads ferromagnetism. For later 
reference we note that Ft is dimensionally an energy 
times a volume which, in this unrenormalized theory, 
is on the order of a Rydberg times a Fermi volume, or 
Ft « eVfci- 



with dij from Eq. (j2.7bl) . The first term has the same 



form as the exchange interaction, Eq. (13.41) . but is much 
smaller, as « Ft(wf/c)^. The second one is the 
dipole-dipole interaction between the electron spins; if 
one replaces the electronic spin density by its quantum 
mechanical and thermal average one recovers the classical 
dipole-dipole term in Eq. ()2.4ap or (|2.7ap . 



3. 0{l/c^): Dzyaloshmsky-Moriya and related interactions 



2. 0{l/c?): Dipole-dipole interaction 

We now turn to terms of 0(l/c^). We first consider 
the vector potential A, which couples to the spin density 
via the B ■ ris term in Eq. Since the coupling is 

directly to the spin channel, no phase space decomposi- 
tion is necessary and integrating out A proceeds as in 
the classical case, except that the A-propagator now is 
frequency dependent, see Eqs. p.2p . We will comment 
on the consequences of this frequency dependence in Sec. 
IIII C 21 below. Neglecting the dynamical aspects of the 
dipole-dipole interaction for now, we obtain a contribu- 
tion to the effective action 



'5'd-d 



27r/^B / dx ns{x) ■ ns{x) 



dxdy S{t^ - Ty)nl[x) dij{x - y) nl{y), 



We now return to the effects of integrating out the 
scalar potential (p. To 0(l/c^) the relevant contribution 
comes from the cross-term that multiplies the Coulomb 
{ipn) term and the spin-orbit (last) term in Eq. (13.11) . 
Contracting the two scalar potentials, integrating by 
parts, using Eq. (jC2b|) . and keeping only terms that are 
bilinear in phase-space spin-density fluctuations, we ob- 
tain a contribution to the effective action 



>S's-o — 



^iim f-ijk / dxdy Vix- y) ijj^^ {x) a^^^ 



(3.6) 



(3.5) After a Fourier transform, Eq. (|3.6I) can be written 




X '^<T3<T2 ) • [(? X p) ^Pa, (k - q/2) (P + 9/2) {p - q/2) (fc + q/2) 



-{qxk) 4,,, {k + q/2) {p - q/2) {p + q/2) ^J,, (k - q/2)] . (3.7) 

Here k = (fc, iujn) comprises a wave vector k and a fermionic Matsubara frequency a;„ = 2T:T{n + 1/2), and p and q 
are used analogously. At this point we generalize to an effective interaction amplitude Vk.pici) that depends on k and 
p in addition to q. Such a structure is generated in perturbation theory from the bare theory, where the interaction 
amplitude is simply given by the gauge field propagator, as we demonstrate in Appendix [d1 We then have 

'S's-o = ^ ( j (9 X p)+VpM{-q) {q x fc)] • {o-^.^^ x cr^3^J 

q.k,p 

x^,, {k - q/2) {p + q/2) {p - q/2) (fc + q/2). (3.8) 

Hermiticity requires pi^) — '^k.p{—q) (see also Eqs. (|B3p ). We now employ the phase space decomposition 
explained in Appendix IB] and focus on the large-angle scattering term, Eq. (jB4bp . Projecting again on the spin 
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density we obtain 



where 



Dig) = y^gkgpik X p) [D(^k+p+q)/2,(k+p-q)/2il«-p)+Vl^k+p-q)/2,(k+p+q}/2ip-k;)] 

T 

^2V 



9k 9p{k-p)xq [Vl^k+p+q)/2,(k+p-q)/2{k - p) - 2'(fc+p-q)/2,(fc+p+q)/2(P ^ k)] 

I 



(3.9a) 



(3.9b) 



Here 



where 



gk = ep <5(efc - ep) 



(3.10) 



is a function that resuhs from the projection onto the 
spin density and pins k and p to the Fermi surf ace. 
Equation (|3.9al) has the form of the resuh obtained by 
Moriya.- In what follows, we discuss the nature of the 
vector D{q) in some more detail. 
£)'^^' has the form 



with d ^'^^ a real vector given by 



(3.11a) 



d^"^^ = iprr'^gkgpip-k) ['D(^k+p)/2,(k+p)/2{k - p) 



2iV 



-^(fe+p)/2,(fc+p)/2(P - k)] . (3.11b) 



D^^) can be written 

Dl'\q) = iD,,,q,+0{q'), 
with Dij a real rank-2 tensor given by 

Dij = Qk 9p <^ilmklPm £j{k,p) , 
k,p 

where 

£i{k,p) = — 



(3.12a) 



(3.12b) 



dxdy {x, - y^)e-'^k+p){x+v)/2 
X j dz sm{{k - p) ■ z)V{x,y;z) (3.12c) 



with ^{x, y; z) the Fourier transform of 'Dk p{q) in anal- 
ogy to Eq. (jB2bp . D^j has a symmetric part and an 
antisymmetric part. The latter can be combined with 
£>'^' above to form a contribution to D that we denote 

by 



D^-\q)=iq xd + 0{q^), 



(3.13a) 



,(2) 



T 
2iV 



Y9k9p{Pikj - k,pj)£j{k,p). (3.13b) 



k,p 



The symmetric part can be written as a diagonal ten- 
sor plus a traceless rank-2 tensor, and the latter can be 
diagonalized by means of an orthogonal transformation 
that amounts to a spatial rotation. The symmetric part 



of Dji we thus can write 



(3.14a) 



(no summation convention) where the obey 

= 0. (3.14b) 

i 

If desirable, the can be explicitly constructed from Eq. 
(j3.12bp . We thus have a second contribution to D that 
we denote by 

ayQy I + 0(q3), (3.15) 



and 



r>(q)-D(+)(g) + £>(-)(g). 



(3.16) 



Combining our results, and transforming back to real 
space, we now have 



where 



12 



Ss-o — iS'dm + 'S'dm + "^DMi (3.17a) 
trZ3 J dx ns{x) ■ (V x ns{x)) , 



(3.17b) 



S'r 



DM 



Qll 

'-'dm 





1 axdx \ 


j dx ns{x) ■ 


Uydy X ns{x) 







(3.17c) 

^ J dx (V-risix)) {d-n,{x)). (3.17d) 
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S'dm is the Dzyaloshinsky-Moriya interaction that is be- 
heved to be responsible for the hehmagnetism observed 
in MnSi and FeGe.— S{y^ is a closely related term but ab- 
sent in systems with cubic lattices (e.g., MnSi or FeGe) 
due to the constraint, Eq. (j3.14bl) . Finally, S^j^ is an- 
other term that is allowed by symmetry, and is generated 
by the above derivation. All three of these terms are con- 
tained in Moriya's general result,^ which takes the form 
of our Eq. (|3.9al) above, but the effects of S^y^ and 
have, to our knowledge, not been discussed explicitly. 

Note that a necessary condition for any of these inter- 
actions to be nonzero is that the system is not invari- 
ant under parity: Both D^^^ and D^^^^ can be nonzero 
only if Vk^piq) is odd under q — !■ q, or, equivalently, if 
'D{x,y;z) is odd under z — > —z. This implies that the 
DM interaction requires a lattice that is not invariant 
under spatial inversion; in any continuum model, where 
the electron-electron interaction is necessarily even un- 
der parity, it vanishes. See Appendix [P] and Sec. II VI for 
further discussions of this point. 



with 

^ijix-y) = S{x -y)Tt+ HbS{tx -Ty)dij{x -y) 



i tr D e^fcj dk + ^ e^kj {ad)k + ^ di dj 



(3.19b) 

Here {ad)k = cikdk (no summation convention). We 
now decouple Eq. p.l9ap by means of a Hubbard- 
Stratonovich transformation with a bosonic field M{x). 
Neglecting a constant contribution to the action, this al- 
lows us to write 



1 



dx dy Mi (x) (a; - y) Mj (y) 



+ - I dxdy [M,{x) T,j (x - y) ni (y) 



-nlix)r.,,{x~y)M,{y)] 



(3.20) 



C. Fermionic action, and magnetic order parameter 

We now have the following result. After integrating 
out the photons, the effective fermionic action reads 



ScsbP, ip] = So + S'i'jjt + Sex + Sd-d + Sum 



(3.18) 



Here 5*0 describes non-interacting electrons (either free 
electrons or band electrons, depending on the model con- 
sidered), and 5;'^^ contains all interactions between modes 
other than the spin density, which we have not explic- 
itly considered with the exception of the Coulomb in- 
teraction, Eqs. (|3.3p . Collectively we denote these two 
terms by £"0. S'ex, <S'd-d, S'dm, Sum, and S^^ are the 
exchange, dipole-dipole, and Dzyaloshinsky-Moriya in- 
teractions given by Eqs. p.4p . (|3.5p . and (I3.17p . respec- 
tively. Sex and Sd-d are always present; which, if any, 
of the terms in S'dm are nonzero depends on the details 
of the lattice structure and absence of spatial inversion 
symmetry is a prerequisite for any of them to be nonzero. 
Collective, we denote the sum of these interactions in the 
spin-density or triplet channel by S'j^^. 



1. Structure of a magnetic order parameter description 

For applications such as fermionic cold gases one will 
want to work directly with the fermionic action. For 
applications to magnets it is convenient to introduce a 
composite field M{x) whose expectation value is propor- 
tional to the magnetization. To this end we write 

Sint dxdy nl{x) F^ (x - y) ni{y), (3.19a) 



If one neglects the interacting part of S'q this action de- 
pends only bilinearly on the fermion fields, and one can 
formally integrate out the fermions in order to obtain a 
theory entire in terms of the order-parameter field M. 
This is a generalization of, and replaces, the Hertz-Millis 
theory.^"'^-'^ However, in general this is not a good strat- 
egy since it amounts to integrating out soft excitations, 
which means that any order parameter theory will in gen- 
eral not be well behaved. Physically, these soft quasi- 
particle excitations can change the nature of the phase 
transition, or they themselves can become criticalj2^ In 
either case they must be treated on equal footing with 
the order parameter fluctuations. It therefore is techni- 
cally advantageous, and physically more transparent, to 
work with the coupled field theory represented by Eq. 
({3:2011 . 



2. Comments on the magnetization dynamics 

We conclude this section with a brief discussion of the 
dynamics of the dipole-dipole interaction, which we ne- 
glected in Sec. IIIIB2I If we restore the frequency depen- 
dence of the photon propagator and expand in powers of 
the frequency, then the leading dynamical contribution 
of the dipole-dipole interaction to Eq. p.20p takes the 
form 



-rdyn 
'd-d 



Mb 



-E 



■ Mi{k)Mj{-k). 



(3.21) 

As long as r2„ scales as |fc|, this scales the same as the 
|r2„|/|fc| term in the order-parameter theory that is in- 
duced by Fermi-liquid effects,-i2, but has a prefactor that 
is smaller by a factor of (wp/c)^. However, in classi- 
cal dipolar magnets the order parameter is known to no 
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longer be conserved.'* That is, scales as a constant, 
and this should be reflected in the quantum theory as 
well, although it is currently not known how this is re- 
alized. This suggests that the contribution shown in Eq. 
p.2ip dominates the Fermi-liquid-induced dynamics in a 
scaling sense, although it has a small prefactor, and will 
become important at sufSciently long time scales. 



IV. DISCUSSION, AND CONCLUSION 

We now discuss the significance of various interactions 
for a number of problems. 



A. Energy scales, and the significance of the 
dipole-dipole interaction 

As we mentioned in Sec. IIII Al the energy scale for 
the exchange interaction in the bare theory is the atomic 
scale, or roughly 100, 000 K. The corresponding length 
scale is on the order of 1 A. This is not consistent with 
the experimental fact that magnetic ordering is observed 
only at much lower temperatures; e.g., on the order of 
1, 000 K in Fe and Ni, on the order of room temperature 
in FeGe, and on the order of 30 K in MnSi. The reason 
for this discrepancy lies in the fact that the bare the- 
ory is renormalized in quantitatively substantial ways, 
and the corresponding energy and length scales in the 
properly renormalized theory are consistent with experi- 
mental observations^ 

Equations (|3.5p and p. 171) show that the dipole-dipole 
and DM interactions, within the framework of the bare 
theory, are weaker than the exchange interaction by a 
factor of (wf/c)^, or about 10~^3^ Relative to the bare 
exchange interaction, this implies an energy scale on the 
order of 10 K, which is comparable with the ordering tem- 
perature in MnSi. On the other hand, the length scale 
associated with the DM interaction (i.e., the pitch length 
of the spin helij*^!) is only about 200 A (in MnSi)^^ to 
700 1 (in FeGe)^= or only a factor of 10^ to 10^ larger 
than the atomic length scale. 

These observations indicate that there are strong 
renormalizations, due to band-structure and many-body 
effects, of all terms in the bare action, and that differ- 
ent terms are renormalized in different ways in differ- 
ent materials. While this makes it hard to make general 
statements, the bare theory suggests that the DM inter- 
actions and the dipole-dipole interaction are generically 
comparable in strength, and in MnSi, for instance, both 
are expected to be a substantial fraction of the (greatly 
reduced by renormalizations) exchange interaction. We 
thus conclude that there is no a priori reason to neglect 
the dipole-dipole interaction in any system where the DM 
interaction is known to be important. This calls for a re- 
evaluation of a number of interesting problems, some of 
which we list in the following subsection. 



B. Significance of the dipole-dipole interaction 

We conclude by discussing a number of problems where 
the dipole-dipole interaction is either known to be im- 
portant, or might be important, with an emphasis on 
low-temperature magnets and other fermion systems. 

(1) Classical Heisenberg ferromagnets. This problem 
was worked on in great detail by Aharony and Fisher— 
for the static critical behavior and by Frey and Schwabl^ 
for the dynamical critical behavior. The renormalization 
group done by Aharony and Fisher started from a non- 
local order parameter theory, Eqs. (|2.7p . which leads to 
a somewhat nonconventional renormalization procedure. 
The nonlocality in Eqs. (|2.7|) is due to integrating out the 
gauge field fluctuations or photons. It would be interest- 
ing to repeat this calculation starting from the coupled 
local field theory given by our Eqs. (|2.2[) before these 
fluctuations are integrated out. 

(2) Classical helimagnets. The standard phase tran- 
sition treatment for helimagnetism due to the DM in- 
teraction is due to Bak and Jensen<^ Neglecting the 
dipole-dipole interaction term they conclude that there 
is a fluctuation-induced first order phase transition from 
paramagnetism to helimagnetism. An interesting ques- 
tion is whether or not the dipole-dipole interaction mod- 
ifies this conclusion. This seems especially relevant for 
MnSi where the phase transition is at low temperatures, 
i.e., it is a weak helimagnet. 

(3) It has been shown that in clean itinerant ferromag- 
nets, the ferromagnetic transition is generically of first or- 
der at zero temperature. ^^'■^'^ This conclusion ignores the 
effects of the dipole-dipole interaction terms. It would be 
very interesting to investigate if dipolar interactions can 
modify this generic conclusion. 

(4) Phase ordering is an important problem in 
ferromagnets.'^^ The dipole-dipole interaction terms has 
not been included in either the classical or quantum 
(zero temperature) ferromagnetic phase ordering prob- 
lems. Simple arguments indicate it will be important. 

(5) Fermionic cold atom systems. Recently there has 
been a considerable amount of work on gases of fermions 
with dipolar interactions, see Ref. jl3i and references 
therein. These systems are important for fermions in 
optical lattices. The dipolar interactions also serve as 
a mechanism for liquid crystal like phase formation in 
fermion systems. 

(6) The dipole-dipole interaction term is important in 
the dynamics of classical antiferromagnets, both in the 
ordered phase, and near or at the phase transition if the 
systems is below its upper critical dimensioui^ Simple 
considerations suggest that they will also be important 
in low-dimensional (1-t-l or 2-1-1) itinerant quantum anti- 
ferromagnets; see Sec. IIII C^ above for one aspect of this 
problem. 
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Appendix A: The complete action to 0{\/(?) 

In this appendix we give the complete action for 
electrons interacting with electromagnetic fields in the 
weakly relativistic limit, up to and including terms of 
0{\/(?). Let — [if,— A) be the 4-vector potential, 
cr = (ci, (72, (T3) the Pauli matrices, and /iB = e/2moC 
the Bohr magneton. We use standard relativistic nota- 
tion, with covariant and contravariant indices related by 
a Minkowski metric 5^1^ — (+,—,—,—). Expanding the 
Dirac equation in powers of 1/c, one obtains, to order 
the following Hamiltonian in first quantization,—"— 

H^Hp+H^o + H^ + Hsm. (Ala) 
Here Hp is the Pauli Hamiltonian, 

Hp = (-iV --A{x,t)Y + eip{x,t) 

Zrric V c / 

-fiBcr ■ {V X A{x,t)) , (Alb) 



and 



^^so — 



ze 



— e 



^cT-{V^{x,t)xV), (Ale) 



V^^ix,t), 



-1 



(Aid) 
(Ale) 



describe the spin-orbit interaction, the Darwin term, and 
the relativistic mass correction, respectively. Via stan- 
dard techniques r^ii^ this theory can be reformulated 
in terms of an action that depends on fermionic (i.e., 
Grassmann-valued) field tJj and its adjoint tp as well as 
the 4- vector-potential field A^. For the partition function 
one obtains 



Z = ZaZ_ 



FP, 



where 



(A2) 
(A3) 



with an action 

SA[^,4';A^]=S4,[^,ij]+SA[Af,] + SS,t,M (A4a) 
Here 

1 _o 



-dr 



2mp 



1 



M^) (A4b) 



describes the electrons with chemical potential and 



Sa[A^] = ^ J dx Af,{. 
, 1 P-1 



dx -dr + V -Aix) (A4c) 
in p J [c J 

describes the electromagnetic fields, with p e 7?, a gauge 

fixing parameter. We use a 4-vector notation x = (t, a;). 



J dx = J dr J dx for space and imaginary time. Sc de- 
scribes the coupling between the fermions and the elec- 
tromagnetic field; it contains four separate contributions: 



Sc[ip,-ip;A^ 



;,L + <JD 



Su + Ss, 



(A4d) 



Here 



Sc,p[i>,ip;A^] = -ie J dx ip{x)n{x) 

+AiB j dx B{x) ■ ns{x) (A4e) 

is the Coulomb and Zeeman paramagnetic coupling 
that is included in the Pauli equation, with n{x) — 

i>a{x)ll)a{x) and ns{x) = IpcrAx) (Ta-j^cr2 "0^2 (2:) • 



Sc,l[iI^,^;A^] = -i2p.B J dx ipcr{x) A{x) -V 
+*A'B J dx CV ■ A{x))n{x) 

dxA^{x)n{x) (A4f) 



2meC^ 

is the Landau diamagnetic coupling; 



Sd[iI^,'4^;A^] = 



dx (V^(/3(x)) n{x) (A4g) 



is the Darwin term that, in the relativistic hydrogen 
atom, leads to the so-called zitterbewegung , and 



dx ipaAx) cr„^c 



— e 
4m^c^ 

■{Vv{x)xV)i;M (A4h) 



is the spin-orbit coupling. 

The second factor in Eq. (jA2|) is 



Zfp = / £) [77,77] e 



'SfP [>):>)] 



(A5a) 



Here 77 is a one-component Grassmannian field known as 
a Fadeev-Popov ghost field, with fj its adjoint, that is 
governed by an action 



'S'fp[77, 77] = dx ri{x)d^d'^ri{x). 



(A5b) 



Appendix B: Phase-space decomposition of 
interaction terms 

For completeness, in this appendix we briefiy recapitu- 
late the arguments that lead to the generation of a spin- 
spin interaction, Eq. (j3.4p . from a density-density inter- 
action, Eqs. (13.31) . For further discussion, see Refs. [22l . 
and 35, 23; the latter also explain the relation to the work 
by Shankarf2S 

Consider an electron-electron interaction with an inter- 
action amplitude W. For simplicity we assume that the 
interaction is purely static, and translationally invariant, 
but otherwise general. The action has the form 
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S ^ ^ J dXi . . . dX4 J dr W{xi - X2,X3 ~ Xi;{x3 + Xi- Xi - X2)/2)Ta,a2a3Cri1pai{xi,T)^„^{x3,T) 

X1pa4{x4,T)^p^^{x2,T), (Bl) 

with r a general rank-4 tensor. We define Fourier transforms 



Wk,p{q) = J dxdydz W{x,y;z), 



(B2a) 
(B2b) 



witli kx — k ■ X — LUnT wlrcre A; is a wave vector and uj„ ~ 2TrT{n + 1/2) a fermionic Matsubara frequency. Hermiticity 
requires 



We tlien have 



k.p.q 



(Ti (T2 (73 <74 



V^ai (fc - q/2) i^a, ijp + g/2) il^a, {p - 9/2) Va. {k + q/2) 



(B3a) 
(B3b) 



(B4a) 



-1 fT 



2 V 



'T X! ^(p-fe+q)/2,(fe-p+q)/2(P + fe) ^^i <T2 <T3 <T4 V'^i ("^ + 9/2) V'<T3 + g/2) 



E VK(fc+p_q)/2,(fc+p+q)/2(p- fc) To-i ^3 ^4 V'.Ti (fc " 9/2) '0<T3 (P + 9/2) V^^i + 9/2) V'ffa (P " 9/2) 



k,p,q 



(B4b) 



2 VI/ 



k,p,q 



xij^,{-p + q/2)^,,{p + q/2). (B4c) 



As long as all wave vectors are summed over, all three 
of these expressions are identical. If one restricts the 
summation in such a way that both |q| and the modulus 
of the third argument of W are smaller than a cutoff wave 
number A, then we can represent the action as a sum of 
all three terms. They represent small-angle scattering, 
large-angle scattering, and 2fcF-scattering, respectively. 
Alternatively, if one is interested in only one of these 
channels, one can pick the appropriate formulation of 5, 
restrict oneself to small wave numbers, and neglect the 
other channels. For our purposes, we are interested in the 
large-angle scattering channel, Eq. (IB4bl) . By choosing 
Ten <T2 era = CT^^ en ^ud making use of Eq. (|C2aP 

we obtain a term that has the structure of the exchange 
interaction, Eq. (|3.4p (in addition to a contribution to the 
number-density interaction). The resulting spin-triplet 
mode is more complicated than a pure spin density, but it 
has an overlap with the spin density and can be restricted 
to the latter by the projection technique explained in Ref. 
I35I . Note that a repulsive Coulomb interaction results in 
an attractive exchange interaction due to a commutation 
of fermion fields that is necessary to write the result in 
the form of Eq. ((Xil) . 

Similarly, by choosing r, — 



making use of Eq. (jC2bp we obtain the structure found in 
Sec. HUB 31 (in addition to terms that couple the number 
density and the spin density). 



Appendix C: Properties of Pauli matrices 

Here we give some properties of the Pauli matrices that 
were used in Sec. IIIIl Let cr = {a^,(jy,a^) = (cr-^, ct^, cr'^) 
be the Pauli matrices with the commutator property 



la'^ (^,J,fc- 1,2,3), 



(CI) 



and cr*^ the 2x2 unit matrix. Then the following identities 
hold: 



1 1 

— _ rr" -I- - /t 



2 
2 



(C2a) 



(C2b) 



' C7 1 (72 (73 (T4 



c^^i <T2 '^'^3 and 



Equation (|C2ap is easily checked by a direct calculation, 
and Eq. ()C2bp follows by multiplying Eq. (jC2aP by cr 
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FIG. 2: Bare interaction V and effective interaction W. 
from the right and using Eq. (|C1|) . 

Appendix D: Nonlocal electron-electron interaction 

In this appendix we demonstrate that a general inter- 
action ampUtude of the type used in Eq. (IBip is gener- 
ated in perturbation theory from an ordinary two-body 
interaction. For simphcity, we consider spinless fermions 
interacting via an interaction potential V; spin depen- 
dence or gradients, as in the last term of Eq. p.ip . are 
easily added. We also consider a translationally invariant 
model; we will comment on this feature below. 

The bare electron-electron interaction is described by 
a term in the action 

I '^^'^y ~ y^'^-'^^^'^^y^'^^y^'^^^^- ^^^^ 

Our goal is to construct an effective interaction of the 
form 

5'off = j dxi . . . dx4 

xW{xi-X2, xz-xa; {xz+xa-xi~X2)/2) 
X'0(xi)'0(x3)V'(a;4) -0(2:2) (D2) 



where G{x — y) = {■ip{x)tp{y)) is the electron propagator. 

In order for W to contribute to Dij we must have 
W{x,y] -z) = —W{x,y;z), see Eq. (|3.12cp . From Eq. 
(ID4p we see that this is the case if and only if V{x) has 
an odd component V~{—x) = —V{x). As far as the con- 
tribution to Dij is concerned we can thus replace V in Eq. 
(|D4I) by V^^, and this automatically ensures W{y,x;z) — 
—W{x,y;z). If V is the propagator of a scalar field, 
such as the quantity V in Sec. IIIIl then these symme- 
try properties can obviously not be realized. However, 
in a realistic solid-state model V represents the screened 



that has a structure necessary for contributing to the 
tensor Aj, Eqs. (I3.12b[ [XT^ . see Fig. H The bare 
interaction corresponds to 



W^'\x,y,;z)=5{x)Siy)V{z). (D3) 



This does not contribute to Dij since it enforces x — y ~ 
0. More generally, contributions to W with the prop- 
erty W{x,y,;z) = W{y,x;z) do not contribute to Dij. 
At second order in V, there are several diagrams, both 
tree-level and one-loop diagrams, that have this property. 
Consider, however, the diagrams shown in Fig. [31 They 




FIG. 3: A second-order contribution to W. 



correspond to 



(D4) 



Coulomb interaction, V{x,y) = vc{x — y)/e{x,y), and 
the dielectric function e has a contribution from the lat- 
tice in addition to an electronic contribution. The latter 
will only have the symmetry of the space group of the 
lattice, and on a lattice without inversion symmetry one 
will have e(y,x) ^ e{x,y). Translational symmetry will 
also be broken, of course, but this is not necessary for 
making the DM interaction nonzero, as the above exam- 
ple shows. Note that one can consider a coarse-grained 
continuum theory for the DM interaction, see Eq. (12.81) . 
but the relevant Landau coefficient must depend on the 



W^'^\x,y;z) = S{x)V{~y) J dx' G{-x' + y/2) G{x' + y/2)V{-x' + z/2) 

+S{y) V{~x) J dx' G{-x' + x/2) G{x' + xl2) V{-x' + z/2), 
I 
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underlying lattice and vanish if the true continuum limit 
is taken. 
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